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Charge carrier localization in extended atomic systems has been described previously as being 
driven by disorder, point defects or distortions of the ionic lattice. Here we show for the first time 
by means of first-principles computations that charge carriers can spontaneously localize due to a 
purely electronic effect in otherwise perfectly ordered structures. Optimally-tuned range-separated 
density functional theory and many-body perturbation calculations within the GW approximation 
reveal that in trans-polyacetylene and polythiophene the hole density localizes on a length scale 
of several nanometers. This is due to exchange-induced translational symmetry breaking of the 
charge density. Ionization potentials, optical absorption peaks, excitonic binding energies and the 
optimally-tuned range parameter itself all become independent of polymer length as it exceeds the 
critical localization scale. Moreover, lattice disorder and the formation of a polaron result from 
the charge localization in contrast to the traditional view that lattice distortions precede charge 
localization. Our results can explain experimental findings that polarons in conjugated polymers 
form instantaneously after exposure to ultrafast light pulses. 


Spatial localization in extended systems has been a 
central topic in physics, since the pioneering work of An¬ 
derson [T] and Mott [2] , and more recently in the context 
of many-body localization [3] . It also forms an important 
theme in materials science of extended conjugated sys¬ 
tems where the dynamics of charges carrier are described 
in terms of localized polarons. unu]. One way to iden¬ 
tify charge localization is through the dependence of its 
energy (e.g., ionization potential or electron affinity) on 
the system size L. In ID systems, if the charge remains 
delocalized, then according to a simple non-interacting 
picture, its energy converges to the bulk limit as l/L^ 
with (a = 1 for a metal or a = 2 otherwise. However, 
if the energy becomes independent of L for L > 4, it 
could be due to charge localization within a critical length 
scale 4- 

Charge localization in conjugated systems can occur 
in several ways: Attachment by point defects [9], lat¬ 
tice disorder effects 13 HO], and formation of self-bound 
charged polarons and neutral solitons by local distortion 
of the nuclear lattice pTUIl] . However, it still remains an 
open question whether localization can occur in disorder- 
free transitionally invariant systems. This question has 
received much attention recently in the context of many- 
body localization p^hUTS] , 

In this letter we provide first-principles computational 
evidence for a new mechanism of localization in ID con¬ 
jugated systems, in which the electrons form their own 
nucleation center without the need to introduce disor¬ 
der into the Hamiltonian. This challenges the widely ac¬ 
cepted picture in which the electronic eigenstates localize 
only after coupling with the lattice distortion [19]. To 


illustrate this mechanism, we study the electronic struc¬ 
ture and the charge distribution in large one-dimensional 
systems with ideal (ordered) geometries. We focus on two 
representative conjugated polymers, trans-polyacetylene 
(tPA) and polythiophene (PT), with increasing lengths 
L = Mil up to M = 70 and M = 20, respectively {ii 
is the length of the repeat unit). Besides their practi¬ 
cal significance |6|, they also exhibit interesting physical 
phenomena, in which polarons, bipolarons and solitons 
affect charge mobility and localization HUaiMll]. 

In Eig.j^we plot the ionization potentials (IPs) for both 
the tPA ( panel a) and PT (panel b) polymers as a func¬ 
tion of the number of repeat units, M. To illustrate the 
effect of localization we focus on the ionization potential 
representing the energy of positive charge carrier (hole) 
rather than on the electron affinity representing the en¬ 
ergy of negative charge carrier (electron), since we find 
the former to localize on shorter length scales (see below). 
Several levels of theory are used: Hartree-Eock (HE) the¬ 
ory, density functional theory (DET) within the local 
density approximation (EDA) [31], the optimally-tuned 
BNL* [321434] range-separated hybrid functional [35|, and 
the B3LYP [36] approximation, and, finally, the GqWo 
many-body perturbation technique (on top of of 
LDA implemented using stochastic DET [38]) within the 
stochastic formulation (sGW) [39]. The LDA and to 
some extent the B3LYP lack sufficient exact exchange 
while HE lacks correlations and screening effects. BNL* 
provide a systematic description of correlations and ex¬ 
act exchange through the process of optimal tuning m- 
GqWo is based on many-body perturbation theory and 
includes exchange, correlation and screening effects m- 
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Figure 1. Ionization potentials (estimated using highest oc¬ 
cupied eigen-energies sh) for (a) trans-polyacetylene and (b) 
polythiophene shown against the inverse number of repeat 
units M in the respective polymer. The repeat unit for each 
polymer is illustrated in the corresponding insets (C, H and S 
are shown by black, white and yellow spheres, respectively). 
Results obtained from different computational approaches are 
indicated by colors and labelled in the figure. Experimental 
data for the ionization potentials (gray circles) were taken 
from Refs. 12311251 and references therein. The dashed lines 
represent a numerical fit to —en (M) = —sh (oc) + ^ for 
LDA and B3LYP (sh (oo) and As are fitting parameters) 

and to —£h (M) = —sh (oo) + As exp M /for HE, 

BNL*, and GW. 


The LDA and B3LYP computations yield IPs that are 
considerably smaller than the experimental values, con¬ 
sistent with previous computational studies on shorter 
polymer chains mnn] and with general theoretical ar¬ 
guments [33 [44] . These IP values drop to their bulk limit 
= 4.4 eV for tPA and = 4.6 eV for PT) asymp¬ 
totically linearly as M~^ for the range of sizes studied 
(they do not fit the purely non-interacting asymptotic 
dependence of In contrast, HE IPs are signif¬ 

icantly closer to the experimental values, deviating by 
less than 0.4 eV. The HE IPs also initially drop as poly¬ 
mer size increases but for a polymer of length exceeding 
a critical value, they quickly converge to an asymptotic 
value, hinting at localization of the hole. The asymptotic 
HE IPs and HE critical length scale are = 6.1 eV 
and = 4.9 nm and = 6.4 eV and = 3.1 nm 
(see Supplementary Material for the approach used to 
determine these quantities). The computational IPs of 
BNL* and the sGW are in even better agreement with 
the available experimental data than those of HE. They 
too show a localization transition with = 5.9 eV 

and = 7.9 nm for tPA and = 6.4 — 6.5 eV 
and = 4.2 — 4.3 nm for PT. Using the results for 
intermediate polymers (which do not exhibit localization 
yet) we can linearly extrapolate to the limit M ^ oc 
and estimate the value of ionization potential if no local¬ 
ization occurs; this yields IP values smaller by 0.5 eV 
which can be viewed as the energy of spontaneous local¬ 
ization. While the asymptotic values of the ionization 
potentials predicted by HE, BNL* and sGW are sim- 
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Figure 2. (a) Calculated optical spectra for selected tPA 

polymers of various lengths (numbers of repeat units M). All 
calculations were performed with the cc-pvTZ basis set using 
TDDFT within the BNL* functional (solid black line with red 
fill) and LDA functional (green filled curve). The fundamental 
band gaps are shown by dashed vertical lines in corresponding 
colors. Red arrows indicate experimental absorption peak 
positions (Refs. [26H29] and references therein), (b) Position of 
the first maxima of the absorption Eopt and the fundamental 
band gap Eg obtained with BNL* and LDA functionals as 
function of inverse number of repeat units. Results for the two 
longest polymers were calculated using the 3-2IG basis set, 
other results are obtained using cc-pvTZ. The exciton binding 
energy is the difference between Eg and the peak maximum 
is illustrated by an arrow. The horizontal full line represents 
the experimental energy of the maximum absorption for the 
infinite system (1.9 eV) [30] . 


ilar, the BNL* and sGW critical length scales £c are 
larger than those predicted by HE. This result is con¬ 
sistent with the tendency of HE to over-localize holes in 
finite systems [45| |46] . 

To further strengthen the validity of the BNL* treat¬ 
ment (and indirectly the GqWo which agrees with the 
BNL*), we compare its predicted optical excitations Eopt 
and fundamental gaps Eg = Sl — £h in tPA to experi¬ 
mental results, where available [26H29] (see Table II of the 
Supplementary material). The absorption spectra shown 
in the left panel of Fig. [^ were calculated using time- 
dependent LDA (ALDA) and BNL* (ABNL*) function¬ 
als [331 ST). It is seen that the ABNL* approach provides 
excellent agreement for the optical gaps in com¬ 

parison to experimental data. The optical gaps Eopt are 
also plotted as a function of 1/M on the right panel of 
Fig. 13 and it is seen that for the largest system studied 
the ABNL* optical gap is in excellent agreement with 
the experimental value [30l [48] . This is in contrast to the 
ALDA results which underestimate this limit by ~ 1 eV, 
and consistently deviate from the ABNL* results as the 
system size increases. In the right panel of Fig. [^ we also 
plot the fundamental gap The values of 

for small systems yields excellent agreement with previ¬ 
ous GqWo results [25] • Furthermore, E^^^* does not 
localize for the tPA lengths studied. Since, local¬ 

izes within a length scale of £c = 7.9 nm the continued 
change in for larger polymers must result from 
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Figure 3. Left panels: The hole densities (top two panels), An(r), for the corresponding labelled methods in long strands 
of M = 50 repeat tPA units (left) and M = 20 repeat PT units (right). The hole is shown as a yellow (aqua) 0.00025aQ ^ 
(-0.00025aQ ^) density isosurface. In the two bottom left panels we plot the cumulative density, p(^), for different functionals. 
The cumulative curve for a shorter tPA polymer with M = 40 (black line) is practically indistinguishable from M = 50. Gray 
areas in the plots show the value of the second cumulant (a) for the corresponding BNL* hole density, which are plotted in the 
right panel for different polymer lengths. The dashed straight line in the right panel is the fully delocalized result (cr = L/\/l2). 
Note that for the larger system we used a smaller basis (3-21G, black symbols) which closely follow the results using a larger 
basis (cc-pvTZ, red symbols). 


a continued change in the electron energy This 

suggests that negative added charge does not yet localize 
for the tPA sizes studied and may explain why the finite 
size gaps are larger than the GqWo gap of 2.1 eV for 
L ^ oo |20l |49l |50] . Note, however, that the GqWo gaps 
are rather sensitive to the size of the unit cell and small 
changes of 0.005 nm in the position of the atoms can lead 
to significant fluctuation of 2.0 to 4.2 eV in the gaps m- 
Since there are no experimental measurements of the fun¬ 
damental gap when L ^ oc, it still remains an open ques¬ 
tion as to the length scale at which electrons localize (as 
opposed to hole localization, which already occurs at the 
system sizes studied). To reach system sizes at which the 
electron localizes will probably require using a stochastic 
approach for BNL* [52]. Finally, panel a of Fig. shows 
that the exciton binding energy Eij = Eg — Eopt is on the 
order of Eg/2 for the larger systems, a value typical of 
other ID conjugated systems [53], indicating that neutral 
excitations are dominated by electron-hole interactions. 

Up to now we have studied localization only from the 
point of view of energy changes. It is instructive to also 
study localization in terms of the hole density, which is 
the difference An (r) = (r) — (r) between the 

ground state density of the neutral (N) and the ground 
state density of the positively charged {N — 1) systems. 
For non-interacting electrons this quantity equals the 
density of the highest occupied eigenstate, which is not 
localized. However, for interacting electrons An (r) must 
be calculated as the difference of densities obtained from 
two separate self-consistent field DFT calculations and 
can thus exhibits a different behavior. We have also as¬ 
certained that the same localization pattern emerges even 
when an infinitesimal charge g ^ 0 is removed, showing 
that localization of the hole density occurs even also in 
the linear response regime. 

The isosurface plots of the hole densities are given in 


the upper left and middle panels of Fig. [^ for the var¬ 
ious methods (excluding sGW). In the lower left and 
middle panels we show the cumulative hole densities 
p (z) = dz' dy^ dx^An (r'). In both types 
of plots it is evident that LDA and B3LYP do not show 
localization of the hole density in any of the systems stud¬ 
ied and in p (z) they show a linear monotonic increase. 
Contrarily, the HF and BNL* charge distributions local¬ 
ize as observed by change of p(z) near the center of the 
chain. In PT this transition in p (z) occurs around one 
of the S atoms closest to the center of the polymer, due 
to the lack of mirror plane symmetry. For long polymers 
exceeding ic^ the BNL* hole density hardly changes as 
seen by the overlapping p {z) of polymers with M = 40 
and M = 50. This implies that the size of the hole is no 
longer influenced by the polymer terminal points and is 
thus independent of system size. 

The extent of hole localization can be described by 
the second cumulant a = yjf An (r') (z' — z)^ dr' (where 
^ = fAn (r') z'dr'). This is shown in the right panel of 
Fig. ^for BNL*. For small sizes a increases as 
consistent with a uniform hole density spread over the 
entire polymer. As L increases beyond £c, the BNL* a 
converge to an asymptotic value, (Too, while those of LDA 
continue to follow the linear L/\^ law (not shown). 

It is important to note that the hole density An (r) is 
dominated by the minority-spin density changes: the or¬ 
bitals having the same spin as the removed electron redis¬ 
tribute such as to localize the hole density near the chain 
center. On the other hand, the majority-spin orbitals re¬ 
main nearly unperturbed and thus do not contribute to 
An (r). This fact reveals that the localization is driven by 
attractive non-local exchange interactions existing solely 
between like-spin electrons. This is further supported by 
the fact that localization only appears in methods that 
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Figure 4. The C-C bond length in the charged M = 50 tPA 
polymer as predicted by HF (left panel) and BNL* (right 
panel) obtained with the 3-21G basis set. In BNL*, a po- 
laron appears as a reduction of the bond-length alternation, 
while in the region about 40 C-C bonds away from the po- 
laron, the alternation is increased to 0.007 nm, similar to the 
experimental value of 0.008 nm for neutral chains |59j . 



account for non-local exchange (HF, BNL*, and GqWo). 

One of the interesting ramifications of the IP stabi¬ 
lization as polymer length exceeds a critical length scale 
is the simultaneous stabilization of the BNL* range- 
separation parameter 7 . This is because in the absence 
of hole localization the tuning criterion, ^ I + Sh = 
0 is expected to become automatically satisfied when 
(semi)local functionals are used in the limit of infinite 
system size [45j I5H456] forcing 7 (and with it the non¬ 
local exchange part of the functional) to drop eventu¬ 
ally to zero. In the systems studied here localization 
saves the day for tuning and the range parameter at¬ 
tains finite asymptotic values of 7 ^^^ = 2.7 nm“^ and 
7 ^^ = 3.1 nm“^. The leveling of 7 with L was reported 
for PT [57], however, it was not previously clear whether 
7 would level-off for tPA. It is worth pointing out that 
7 does not change significantly (< 0.2 nm“^) when LYP 
correlation is used instead of LDA in the BNL* calcula¬ 
tion. 

While HF supports partial localization, its hole density 
also exhibits oscillations along the entire polymer length 
that do not diminish with system size. These indicate a 
rigid shift of charge between neighboring atoms: From 
double to single C-C bonds in tPA and from S to nearby 
C atoms for PT. This is consistent with the tendency of 
HF to eliminate bond-length alternation in the entire tPA 
polymer chain [58] as shown in the left panel of Fig. 
BNL* on the other hand eliminates the bond-length al¬ 
ternation only in proximity of the localized hole density 
(right panel of Fig.|^, consistent with a localized polaron 
model. 

In summary, using first principles density functional 
theory and many-body perturbation theory, we have 
shown that positive charge carriers can localize in ID 
conjugated polymers due to a spontaneous, purely elec¬ 
tronic symmetry breaking transition. In this case, local¬ 
ization is driven by non-local exchange interactions and 
thus cannot occur when (semi)local density functional 
approximations are used. HF theory, which has non-local 
exchange, shows a localization transition in a relatively 


small length-scale but predicts complete annihilation of 
bond-length alternation upon ionization, irrespective of 
polymer length. BNL*, which through tuning includes 
a balanced account of local and non-local exchange ef¬ 
fects, provides an accurate description of the optical gap 
in comparison to experiments and shows a localization 
transition with a length scale (estimated from the lev¬ 
eling off of the IPs) that agrees well with the 5 GW ap¬ 
proach. Moreover, BNL* predicts a localized disruption 
of the bond-length alternation. 

The localization phenomenon is driven by the same- 
spin attractive non-local exchange interactions and there¬ 
fore, cannot be explained in terms of classical electro¬ 
statics. There is no reason to assume that the observed 
emergence of the localization length £c in finite systems 
will not readily occur also in infinite systems, where hole 
states near the top of the valence band are necessarily 
infinitely degenerate. 
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